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SUMMARY 

Several types of propulsion systems are evaluated for the 0. 1-AU solar probe mis- 
sion commencing from a low, circular Earth orbit. The four propulsion systems com- 
pared are (1) the chemical rocket, (2) the nuclear rocket, (3) the electric rocket, and 
(4) the chemical-electric hybrid system. The chemical-electric hybrid system consists 
of a chemical first stage fired until at least escape velocity is attained and an electri- 
cally propelled second stage. For this comparison, calculus of variations calculations 
were performed upon a mathematical model set up for electric and chemical-electric 
hybrid systems employing constant-thrust engines. 

It is shown that the chemical- electric hybrid is by far the most attractive system 
for present and near future powerplant technology. In particular, the Saturn IB booster 
cannot deliver any payload to 0. 1 AU if existing (or near future) hardware is utilized 
regardless of the upper stage propulsion choice except for the chemical-electric hybrid 
system. The study considers mission times between 60 and 650 days and specific 
electric powerplant masses from 10 to 150 pounds per kilowatt. 


INTRODUCTION 

Unquestionably mankind' s environment is predominately controlled by the Sun. It is 
only natural then to increase our scientific knowledge of the "center of attraction" in 
order to understand and anticipate the behavior of our environment and its underlying 
mechanisms. At present, our knowledge of the Sun is based almost entirely upon limited 
experiments performed on Earth. The protective atmospheric and electromagnetic 
shields surrounding the Earth limit Earth-based solar observations to optical and radio 
frequencies in the electromagnetic spectrum and severely restrict the study of solar 
particle emission. 
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This situation is being alleviated in some degree by orbiting solar observatories and 
interplanetary probes. There remains, however, a great deal of desirable knowledge 
that requires a direct close-in study of the Sun (refs. 1 and 2). Requiring a solar probe 
to attain a perihelion radius of 0. 1 AU is justified on the grounds that Sun-probe dis- 
tances of this magnitude are required by some of the important proposed solar experi- 
ments (corona radar sounding, ultraviolet and X-ray spatial resolution of the solar disk, 
determining the connection between the rotating solar corona and the solar wind, etc. ). 
The intent of this note is not to discuss these experiments or the makeup of the scientific 
payloads that perform them, but rather to analyze vehicle systems capable of accomplish- 
ing this mission - with emphasis on the combined chemical- electric system. 

There have been several solar flyby probe studies in the past. The case of high 
thrust transfers utilizing one or more heliocentric pulses is analyzed in references 
3 to 6. The major conclusion from these studies is that only very small payload 
fractions exist (if any at all) for close solar flybys using chemical propulsion. 

An evaluation of the solid- core nuclear rocket for this mission was done in refer- 
ence 6 for Saturn IB class boosters. Although some payload capability was reported in 
the study, it must be considered quite optimistic since the tungsten- core nuclear engine 
was assumed to weigh only about 1400 pounds. (Nerva class engines are about an order 
of magnitude heavier.) The 5000-pound graphite- core nuclear engine assumed in refer- 
ence 7 for this mission allowed no payload for Saturn IB class boosters. 

The application of a constant-thrust electric propulsion system to this mission has 
been studied in some degree in references 6 to 8. In all cases, the specific mass of the 
electric propulsion system was assumed to lie far below present and near future capa- 
bilities. Furthermore, the electric stage was assumed to commence either in Earth 
orbit or at Earth escape, thus avoiding the question of when should chemical propulsion 
cease and electric propulsion begin. 

The present study evaluates constant-thrust electric propulsion in light of both 
future and present capabilities with regard to powerplant specific mass (values between 
10 and 150 Ib/kW are assumed). It also investigates the problem of optimum stagii^ 
between high- and low-thrust stages. In particular, it is assumed that electric pro- 
pulsion begins either in Earth orbit (designated herein as the electric system) or at some 
arbitrary velocity following Earth escape (designated herein as the chemical-electric 
hybrid system). 

A comparison between the various systems is then made in terms of payload ratio, 
mission time, and electric power requirements. The power requirement of the payload 
is an important parameter for this mission since payload ratios are small and the power 
required for onboard experiments and telemetry in the vicinity of the Sun can be quite 
large. Power estimates range from several hundred watts to several kilowatts (ref. 8). 
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ANALYSIS 


I 


No analysis is presented herein for chemical and nuclear systems. Instead, the 
results of other studies such as reference 7 are used where needed. This analysis is 
concerned with determining the performance of electric and chemical-electric hybrid 
systems initially in a 100-nautical-mile circular Earth orbit. Payload ratio will be the 
primary performance criterion. Actually, mission time and electric power are of im- 
portance, too, and an example of the tradeoffs among these three entities will be given 
in the section RESULTS AND DISCUSSION. The chemical-electric hybrid system is 
treated first. 


Hybrid System 

The problem is to determine the maximum payload ratio achievable with a chemical 
first stage and a constant- thrust electric second stage, holding mission time, specific 
powerplant mass, and other vehicle parameters constant. 

To solve this problem with precise results, it is necessary to consider the three- 
body, variational problem with optimum staging. In order to avoid the complex calcula- 
tions involved with such a problem, three assumptions were made: (1) the chemical 
escape assumption - the chemical stage is fired until at least escape velocity is at- 
tained, (2) the two-body assumption - the gravitational effect of the Sun is negligible 
during the chemical stage flight, while the gravitational effect of the Earth is negligible 
during the electric stage flight, and (3) the zero time assumption - the time elapsed be- 
tween chemical stage burnout and sphere of influence penetration is negligible so that 
electric propulsion commences at the Earth’s sphere of influence in heliocentric space. 
The two-body assumption cannot be realistically made unless the chemical escape 
assumption already holds. Although some accuracy is lost due to the two-body and 
zero time assumptions, the main drawback is that chemical stage burnout velocities 
less than escape velocity are not allowed (chemical escape assumption). Nevertheless, 
it is felt that adequate results should be obtained, particularly, when high burnout veloc- 
ities are found to be optimum. 

The thrust vector control of the heliocentric phase is determined by variational 
principles. This means that (1) the Euler- Lagrange equations are employed for the 
determination of the electric thrustor orientation, (2) the transversality relations are 
used to optimize the heliocentric travel angle and Earth escape velocity orientation, 
and (3) the coast phases are optimized. 

The payload equation for the chemical- electric hybrid system just described may 
be written 
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where is the payload ratio of the chemical stage and the sum of terms enclosed 

in parentheses represents the payload ratio of the electric stage. (All symbols are de- 
fined in the appendix. ) In order to facilitate the discussion concerning the maximization 
of the payload ratio mj^M^ equation (1) will be rewritten in terms of the pertinent 
problem variables (specific impulse, initial acceleration, etc.). 

The tankage mass of the electric stage m^ is usually taken to be proportional to the 
propellant mass of the electric stage m . Thus, if k is the proportionality constant, 

Jr 

m^. = kmp (2) 

Defining the specific powerplant mass oi to be the ratio of the powerplant mass 
divided by the power supplied to the thrustors P and defining the thrustor efficiency tj 
to be the ratio of the propulsive power divided by the total power supplied to the thrustors 
result in the powerplant mass fraction being rewritten as 


”^pp _ Q, _ ^o^e^o q_ 
”^o "^0 47.47X10^ 


( 3 ) 


where a is the initial acceleration of the electric stage, I is the specific impulse of 
O q 6 

the electric stage, and 47. 47x10’^ is a constant that is required for the system of units 
employed in this report. 

The payload mass of the chemical stage (initial mass of the electric stage) can be 
written as 


m^ = - M. 

o o p h 

where the hardware mass Mj^ is composed of such things as tankage mass, engine 
mass, guidance and control mass, and structure mass. For chemical rockets, the 
hardware mass can usually be taken to be proportional to the propellant mass M . If K 
is the hardware proportionality constant, then 

= KMp (4) 

Assuming that the chemical stage imparts an impulsive velocity change from circular 
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orbit velocity to burnout velocity the chemical stage payload ratio is given by 


— = (1 + K)e ° ^ 

Mo 


( 5 ) 


If equations (2) to (5) are substituted into equation (1) along with the final mass rela- 
tion m. = m - m„, there results 
lop 


m. 


M, 


(1 + K)e ^ ° - K 




(l + k)-i 




o e°o 


®- (k+ ^ 


m. 


“o 47. 47X10^ 


m. 


( 6 ) 


This is the final form of the function to be maximized. The constants in this equation 

must be ascribed values and these are given in table I. The electric stage final mass 

ratio mr/m increases as a.. I , and Vk are increased. However, increases in 
r o o^ e D 

a^ and 1^ also increase mp^/m^, while increasing decreases m^M^. Clearly, 

there exist values of a^, 1^, and Vy^ that will result in a maximum payload ratio. 

These three variables may be optimized by a three-dimensional search scheme. 

The maximization of m^/m^ can be formulated as a Mayer problem in the calculus 

of variations wherein mission time is treated as a parameter and the constraints are 

composed of the two-dimensional equations of motion plus the constraint requiring the 

thrust magnitude to be either zerp or some constant value. This problem has been 

treated often in the past and its solution is given in reference 9. In the problem at hand, 

the heliocentric travel angle is left unspecified (free for optimization). Another boundary 

condition that is not completely specified is the initial heliocentric velocity V^, which 

is obtained by vectorially adding the velocity 

relative to the Earth at the sphere of influence 

V_ to the Earth’ s orbit velocity about the Sun 
s 

as shown in sketches (a) and (b). The 
optimum orientation of V_ is in the direction 

o 

of the initial electric thrust vector as can be 
shown by the transversality condition of the 
Mayer formulation. The magnitude of V_ 

o 

is related to the burnout velocity Vb by the 
conservation of energy as follows: 


TABLE I. - ASSUMED CONSTANTS 


Parameter 


Assumed value 


Circular Earth- orbit I 25 600 

velocity, V^, ft/ sec 

Specific impulse of I 420 

chemical stage, 1^,, sec 
Hardware fraction of I 0. 137 

chemical stage, K 

Ratio of tankage to propellant | 0. 02 

mass of electric stage, k 
Ratio of structure to initial I 0. 04 

mass of electric stage. 



(7) 
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(a) Earth -centered escape maneuver. 




Vs 


Earth 

(b) Heliocentric initiai velocity Vq. 

For this study, the sphere-of-influence radius r was taken to be 120 Earth radii. 

s 

Finally, the specific powerplant mass a and thrustor efficiency 77 must be con- 
sidered. Payload ratio is a monotonic function of a; hence, a can be treated as a 
parameter. Dealing with the thrustor efficiency is a more difficult problem because it 
is a function of specific impulse, and the functional relation is dependent upon the type 
of thrustor employed and its state of development. One approach to this dilemma is to 
assume some characteristic function 77(1^) applicable to a particular thrustor. This 
obviously leads to results that are dependent on the thrustor choice. Alternatively, the 
choice of thrustor may be delayed by simply treating the ratio a/ri as a parameter in- 
stead of 0! (as in ref. 6). However, this approach ignores the dependence of 77 on 1 ^, 
and therefore does not arrive at the true maximum payload ratio except for a hypothet- 
ical, constant-efficiency thrustor. Both of these approaches were tried and are exam- 
ined in more detail in the RESULTS AND DISCUSSION section. In either case, the pay- 
load ratio may be maximized by a three-dimensional search on Vb’ and I^ over 
the class of optimal heliocentric trajectories for any given pair of mission time and 
a (or a/rj). 



Optimization Of V|j 

Consider a hypothetical case wherein mission time and a have been chosen, both 
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and 1^ have already been optimized, and 
only remains to be optimized. Since 
is constrained to take on values at least as 
large as the escape velocity it might be 
expected that the optimum could fall ex- 
actly on this limit (i. e. , = Vg). How- 

ever, the mathematical model that has been 
created for the hybrid system ensures that 
^ ^ ^ Qpj > Vg. That this is so can be easily 

° '^s,opt s demonstrated with the aid of sketch (c). In 

(c) Optimization of V5 for hybrid system. this sketch typical curves of chemical stage 

payload ratio (m^/Mg), electric stage payload 
ratio (mj^/m^), and overall payload ratio (mj^/M^) are displayed as functions of the 
vehicle velocity at the sphere of influence, Vg. For small values of r^/r^, the vehicle 
is considered to have exactly escape energy when V_ = 0; hence, from equation (7), the 
burnout velocity is exactly escape velocity when 


i|— 



V 


b 




( 8 ) 


Since m^^/Mg = (m the slope of the overall payload ratio curve is 


d 





(9) 


From equations (5) and (7), it is easily shown that the slope of the chemical stage 
payload ratio curve is 


/m \ 

9 |— 1 

JjV - .LtK^Sg-^VVcVlc^o 
Ic^o Vb 


( 10 ) 


Thus, when V„ = 0, the slope of the chemical stage payload curve vanishes and the 
s 

slope of the overall payload ratio curve as given by equation (9) becomes 
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( 11 ) 



The greater than zero inequality can be written since dim^ym^/dV ^ and are 

positive for all cases of interest. Hence, when Vg = 0 (i. ©• > = V^), the overall 

payload ratio curve has a positive slope. 

On the other hand, if Vg is large enough, the slope of the electric stage payload 
ratio curve 9(mj^/m^)/9Vg is bound to be negative. For such a Vg both terms of 
equation (9) are negative and hence the overall payload ratio curve has a negative slope. 

If V^ denotes the value of Vg where 9(mj^/m^)/ 9Vg = 0, Vg has been shown to be 
bounded by. 


0<Vs,opt<V 


and hence an optimum Vy^ exists for which 


Vopt>Ve 


( 12 ) 


(13) 


Electric System 

The payload maximization scheme for pure electric propulsion trajectories may 
be performed similarly to the chemical- electric hybrid system scheme. The burnout 
velocity Vj^ is not involved so that only two variables (a^ and 1^) need to be optimized 
by the search process. Optimum Earth-escape trajectories for low-thrust vehicles are 
very nearly tangential thrust spirals that can be quite accurately computed in closed 
form (e. g. , ref. 10). Thus, only the heliocentric portion of the flight needs to be 
optimized. In this case, the ratio m^M^ in equation (1) is unity, and the Earth- escape 
spiral propellant mass mpg may be computed after the heliocentric optimization from 
the following equation (derived in ref. 10): 

"^ps = ® 


where 
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u = 


(15) 


I g 

e^o 

The empirical correction factor y is taken to be a function of the ratio of a^ to the 

local gravitational acceleration g and is displayed on figure 1 (reproduced from ref. 10). 

The time Tj^ required to spiral to escape velocity is computed from the propellant mass 

and propellant flow rate m__ as follows: 

ps 

T, = v(l - e-‘')^2?2 (16) 

“ps “o 

The aforementioned treatment of the electric system optimization implicitly assumes 
that the initial heliocentric position and velocity are identical with the Earth's heliocen- 
tric radius and orbital velocity. 


RESULTS AND DISCUSSION 
All-Propulsion Constraint 

Preliminary data showed that the optimum payload trajectories tended to be nearly 
all propulsion (only a short terminal coast phase). Furthermore, as will be shown in 
the section Companion of thrustor types , trajectories constrained to all propulsion 
operation resulted in payload ratios that very nearly equalled the payload ratios obtained 
when coast phases were allowed. Since the all-propulsion constraint leads to significant 
simplifications of the problem, it was imposed whenever it was appropriate to do so. 
Two simplifications arise from this constraint: (1) the boundary-value problem associ- 
ated with the variational problem becomes much less sensitive and therefore the itera- 
tive method of solution converges faster, and (2) a^ is no longer a variable that needs 
to be optimized; instead, it is determined by the associated boundary-value problem. 


Payload Ratio 

By using the approach described in the ANALYSIS section wherein a/rj is treated 
as a parameter, and with the all-propulsion constraint imposed, payload ratios were 
calculated and are plotted as a function of mission time and qj/tj in figures 2(a) (hybrid 
system) and (b) (electric system). Lines of constant optimum specific impulse are also 
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shown and are seen to be fairly flat over most of the range of mission times. Moreover, 
the optimum specific impulses are relatively low (2000 to 4000 sec) for the region of 
greatest immediate interest - high a/ri. The data for the chemical and nuclear systems 
appearing in figure 2 were taken from reference 7, wherein the payload ratio was taken to 
be an increasing function of initial spacecraft mass. Hence a small spectrum of payload 
ratios is represented in figure 2 as bars. The tops of these bars correspond to space- 
crafts that could be inserted into earth orbit by Saturn V class boosters. Spacecrafts 
inserted into earth orbit by Saturn IB class boosters cannot deliver a positive payload 
ratio and are not represented in figure 2. The payload performance of chemical and 
nuclear systems is quite unattractive. 

To compare the hybrid and electric systems, figure 2 was superimposed as shown 
in figure 3. For each value of a/ri there is a curve for the electric system and another 
curve for the hybrid systems. These curves intersect at the point where both systems 
deliver the same amount of payload. The solid curve in figure 3 passes through these 
intersection points and thus divides the figure into two regions - one in which the electric 
system offers greater payloads and another in which the hybrid system offers the greater 
payloads. Typical present and near future powerplant systems fall into the region where- 
in the hybrid vehicles offer greater payload capability. 

The existence of a region in figure 3 where hybrid systems are poorer than all- 
electric systems is a consequence of the ^ restriction. As defined earlier, the 
chemical and electric systems are simply limiting cases of the unrestricted chemical- 
electric hybrid system. That is, the hybrid system degenerates into the electric sys- 
tem as Vjj approaches whereas it degenerates into a single-stage chemical system 
if Vjj becomes large enough. Therefore, a plot of payload ratio as a function of chemi- 
cal stage burnout velocity that incorporates a proper optimization (arbitrary in- 
clusion of n-body effects, and inclusion of electric propulsion inside the Earth's sphere 
of influence) might look something like the solid line in sketch (d). The optimized 
hybrid system (unrestricted) can never yield a payload ratio lower than either the chemi- 
cal or electric systems since it includes these systems as degenerate cases. However, 
if the optimization scheme relies on the assumptions of two-body motion and zero electric 
propulsion within the Earth’ s sphere of influence and restricts to values greater 
than Earth- escape velocity, then according to equation (13), an optimum must also 
occur at some value greater than Earth-escape velocity. This is illustrated by the 
dashed line in sketch (d). Of course the true optimum Vj^ can be as low as V^. It is 
only when the true optimum lies in the region :< (as depicted in 

sketch (d) that the hybrid system as treated herein is credited with a payload ratio 
significantly lower than it ideally could deliver. 

One further point, the dashed line representing the approximate solution in 
sketch (d) suggests that the exact solution may be shaped quite differently than its 
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representation in this sketch. Indeed, 
the exact solution may contain an in- 
flection point and perhaps possess two 
optimums. This question remains un- 
resolved at present, but it is suspected 
that the inclusion of electric propulsion 
within the Earth' s sphere of influence 
could remove the con- 

dition and produce a curve of the shape 
depicted in sketch (d) for the exact 
solution. 

To compare the electric and 
hybrid systems for a specific thrustor, 
it is necessary that the thrustor ef- 
ficiency be known as a function of specific impulse. Thrustor efficiency data (taken from 
refs. 11 to 17) is presented in figure 4 for several types of thmstors. As an example, 
the efficiency function characteristic of small electron-bombardment thrustors with 
mercury as the propellant was used in conjunction with figure 2 to produce the payload 
curves shown in figure 5. The sharp changes in the slope of the curves on figure 5 
occur because the lower part of each curve represents a restricted (Vj^ > Vg) hybrid 
system whereas the upper part represents an electric system. That is, for each a 
only the greatest payload ratio was plotted - regardless of which type of system a partic- 
ular point belongs. Presumably, if completely arbitrary were allowed for the 
hybrid system, the payload curves would not suffer a discontinuous slope. 

Comparison of optimization schemes. - The payload capability of small electron- 
bombardment thrustors was also computed using the alternate approach discussed in the 
ANALYSIS section, that is, treating o! as a parameter and including the efficiency 
curve in the optimization scheme. These results are given in figure 6. The data of 
figure 5 are reproduced on figure 6 for a direct comparison between the two approaches. 
The simpler approach, which treats a/t] as a parameter, produces payload curves 
that have the same general shape as the more exact approach. However, in the region 
of prime interest for the near future - high o; and low payload ratio - the error pro- 
duced by the parametric a/r] approach is relatively large. This is because the optimum 
specific impulses for this region lie on the lower end of the r\ against 1^ curves - 
where ?] is a strong function of 1^. Hence, the remaining figures in this report are 
based on data that include the efficiency curve in the optimization scheme. 

Comparison of thrustor types. - In addition to the small electron-bombardment 
thrustor, two other electric engines were investigated; namely, a large (50- cm diam- 
eter) electron-bombardment thrustor and an electrothermal thrustor. The character- 



11 



istic efficiency curves for these engines are shown in figure 4. In the case of the 
electrothermal engine, the tankage fraction was increased to 15 percent since hydrogen 
is the usual propellant for such engines. The results are shown in figure 7. The gain 
in efficiency due to increased thrustor size is reflected by a substantial performance 
increase. The performance of the electrothermal engine is inferior to that of the 
electron-bombardment engine because of the sizable increase in tankage fraction and un- 
favorable efficiency curve. 

Several cases of optimal coast phase trajectories are also shown in figure 7 for the 
large electron-bombardment thrustors. It is evident that very little payload benefit 
results from coasting trajectories - thus justifying the selection of all -propuls ion tra- 
jectories throughout much of the report. 

The optimum specific impulse, high-thrust burnout velocity, and initial power-mass 
ratio associated with the small size electron-bombardment and electrothermal engines 
are presented in figures 8 to 10. As a general rule, it can be stated that, as the mission 
difficulty increases (shorter mission time and/or increased specific powerplant mass), 
the optimum specific impulse decreases while the optimum high- thrust burnout velocity 
increases. Figure 10 shows that the hybrid system power requirement is about one- 
third of the electric system power requirement. 


Payload, Power, Mission Time Tradeoffs 

Certainly payload is an important parameter to be considered when planning a mis- 
sion. Other parameters may be of equal importance, however, so that a realistic 
mission design usually involves a compromise between payload and other factors. One 
such important factor for the 0. 1-AU solar probe mission is the payload power require- 
ment. Nonelectric propulsion system studies must include the mass of a nonpropulsive 
electric power supply (for telemetry, experiments, etc. ) as part of the payload. The 
’’useful payload” is the result of subtracting the power supply mass from the ’’total 
payload. ” The useful payload of electric propulsion systems may be calculated in the 
same way; however, nonpropulsive power often need not be carried. For instance, 
during the terminal coast phase associated with optimal flyby trajectories, the power 
normally available to the electric engines may be available to increase the power avail- 
able for data transmission and payload experiments. Furthermore, in some cases 
where propulsive power and nonpropulsive power are needed simultaneously, the non- 
propulsive power would cost less than its equivalent for a nonelectric propulsive system. 
This is because the specific power supply mass of some systems decreases as power is 
increased and propulsive power is usually many times greater than nonpropulsive power. 

In order to illustrate the tradeoff relation between useful payload, power, and 
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II 


TABLE n. - SYSTEM PERFORMANCE UTILIZING A SATURN IB BOOSTER®^ 


System 

Mission 

Maximum 

Electric 

Specific 

First- stage 


time, 

useful 

power. 

impulse, 

burnout 


’’^M’ 

payload, 

kW 

sec 

velocity 


days 

lb 



ft/sec 

Chemical rocket (Centaur + 

75 

0 

0 

444, 444 


”kick’* stage) 






Nuclear Rocket (Nerva class) 

75 

0 

0 

800, 800 


Electric rocket^ 

400 to 650 

0 







Chemical- ele ctrothermal 

500 

850 

15 

420, 1500 

40 000 

hybrid^ 






Chemical - large electron- 

400 

1020 

35 

420,3200 

41 000 

bombardment hybrid*^ 

500 

2360 

37 

420, 3700 

38 700 


^Saturn IB assumed to inject 32 000-pound payload into 100-n. -mi, parking orbit. 
*^Specific powerplant mass a = 100 Ib/kW. 


mission time, these criteria are shown in figure 11 for typical near future systems. All 
systems were assumed to have been boosted into earth orbit by a Saturn IB and have a 
specific power supply mass of 100 pounds per kilowatt. The two-stage high performance 
chemical system data and the one and one-half - stage nuclear system data is presented 
only at the 75-day mission time since this time represents a minimum energy transfer 
(i. e. , higher mission times are not profitable). Neither of these two systems can 
deliver positive useful payloads - the negative payloads shown in figure 11 are based on 
data from an unpublished study by the Advanced Development and Evaluation Division at 
the Lewis Research Center and reference 7. The pure electric system is also unable to 
deliver positive payloads and is not represented in figure 11. Thus figure 11 shows that 
there is no competition among the various systems for this particular case - only the 
chemical- electric hybrid system can deliver any useful payload. For some other, less 
difficult missions, there will exist a definite competition between the various systems 
and a figure similar to figure 11 would display the tradeoffs between useful payload, 
power, and mission time. In general, high power levels favor the electric and chemical- 
electric hybrid systems. The important point is that the payload power requirement can 
strongly influence a system choice in terms of useful payload. 

Table II summarizes the results for a Saturn IB booster. It is true that if a larger, 
more efficient spacecraft were assumed (which covild be boosted into Earth orbit by a 
Saturn V, for example), chemical propulsion alone could deliver some payload (ref. 7) - 
but by the same token the hybrid system’s payloads would increase correspondingly. 

The nuclear rocket system is similar to the chemical system in that larger boosters are 
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necessary to provide positive payloads. In addition, 
nuclear rocket systems do not fall into the near 
future category. Clearly, if both Saturn V class 
boosters and nuclear rockets were available, the 
nuclear-electric hybrid (instead of the chemical- 
electric hybrid) should be compared to the nuclear 
system. 

Figure 11 and table II illustrate the fact that, 
while neither the chemical system nor the electric 
system can accomplish the mission, the hybrid sys- 
tem can. In fact, specific power supply masses on the order of 100 pounds per kilowatt 
can be realized almost immediately with 5- to 30-kilowatt class solar cell modules 
(refs. 18 and 19) and 400- to 500-day thrustor lifetimes are reported to be within the 
present state of the art (ref. 11). 

The foregoing analysis has been carried out under the tacit assumption of constant- 
propulsion power. However, in the case of solar energy dependent power supplies such 
as solar cell, thermionic, and thermoelectric systems, the power will vary as the 
distance to the Sun changes unless special means are provided to keep the power constant. 
This is due to the combined effect of (1) the solar energy flux varying as the inverse 
square of radius and (2) the power conversion efficiency varying as some function of 
temperature. Hence, the power profile (power as a function of distance from the Sun, 
see sketch (e)) of such systems is a complex function involving the type of system and the 
particular system design as well as the distance from the Sun. In terms of power pro- 
files and specific mass, solar cell power supplies appear more attractive than both 
thermoelectric and thermionic systems. 

One drawback of solar cells is that such a system designed for operation at Earth’s 
radius or greater would experience a sharp power decay between 0. 2 and 0. 3 AU. While 
this is not important from a propulsion standpoint (optimum coast phases begin at 
greater radii), it is quite unfortunate from a payload and communications power stand- 
point. Thus, a small thermionic generator or thermoelectric panel, whose power out- 
put is quite attractive if designed for 0. 1-AU operation, might be carried along to sup- 
ply power to the payload and communication system when the solar cell power output 
vanishes. Of course, the mass of such a nonpropulsive power supply should be sub- 
tracted from the payload to yield useful payload as discussed in the beginning of this 
section. 

The effect of a nonconstant power profile on payload ratio was approximated by com- 
puting an average power over the trajectory and using a payload against power plot such 
as figure 11. K a typical solar cell power profile designed for 1. 0-AU operation were 
assumed, such calculations revealed that, for a 500-day mission employing a 100-pound- 



Distance from Sun, AU 
(e) Typical power profile for solar cells. 
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per-kilowatt power supply (a based on 1. 0-AU exposure), the payload decrease would 
be about 20 percent whereas for a 400-day mission with the same power supply the de- 
crease would be about 10 percent. For shorter missions or smaller values of a, the 
decrease would be less. The existence of an optimum mission time for radius dependent 
power supply systems (suggested by the previous statements) remains open to question 
since variational solutions incorporating variable power would modify the constant power 
trajectories in such a manner that the decrease in payload ratio would be minimized. 


Effect of Perihelion Radius 

Although one of the basic ground rules for this study is the 0. 1-AU perihelion radius 
requirement, it is of interest to see how changing the perihelion radius affects the solar 
probe mission. Figure 12 displays the variation of payload ratio as a function of peri- 
helion radius for a hybrid system employing large electron-bombardment thrusters and 
a 100-pound-per-kilowatt power supply. It is evident that the perihelion radius has a 
strong effect on payload ratio. However, for mission times on the order of 500 days, 
perihelion radii of even less than 0. 1 AU are readily attainable with the hybrid system. 


Trajectory Characteristics 

Effect of mission time. - Two typical trajectories of the hybrid system are shown 
in figure 13. The long mission time trajectories swing far outside Earth's orbit - in 
fact, for mission times greater than 350 days, a Mars swingby maneuver may be attrac- 
tive. Arrows denoting thrust direction are placed on the trajectory at equal time incre- 
ments. It is evident from the trajectory diagram that the long missions are character- 
ized by an initial phase that moves the vehicle into a state of low velocity and a terminal 
phase which simply removes angular momentum. Figure 13 also shows that long mis- 
sion time probes have the potential of gathering data over a broader spectrum of solar 
distances which could at least partially balance the disadvantages of long mission times. 

Terminal configuration angle. - An important parameter in communications studies 
of solar probes is the terminal configuration angle between the Sun-Earth line and the 
Sun-probe line as shown in figure 14. The Sun is a powerfvil source of background noise, 
and there is a communications blackout region when the configuration angle is near 
0° or 180*^. The terminal configuration angle for the missions employing the small 
electron-bombardment engines is shown in figure 14. 

Absolute and local optimal trajectories. - The bulk of the data presented in this re- 
port depends upon the generation of absolute optimum heliocentric trajectories. Such 
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trajectory calculations can be very difficult since they require solutions to boundary- 
value problems. Another pitfall is the occasional occurrence of ’’local optimum” solu- 
tions. This is a consequence of the Euler- Lagrange equations being necessary but not 
sufficient conditions for absolute optimality. During the data generation process, 
optimum trajectories were encountered in the region of 250- to 350-day missions, which 
proved to be local optimums rather than absolute optimums. Two optimum trajectories, 
one local and one absolute, are shown in figure 15 for the same mission. Note that the 
locally optimum trajectory is characterized by a change from posigrade to retrograde 
motion. In this study, trajectory solutions that were not absolute optimums were easily 
identified by this characteristic posigrade-to- retrograde motion and the very poor con- 
vergence associated with their boundary-value problems. 

Heliocentric AV. - In order to facilitate payload computations under weight assump- 
tions different from those assumed in this report, heliocentric characteristic velocity 
AV data is presented in figure 16. Since all- propulsion trajectories are nearly optimum 
and since AV is quite insensitive to specific impulse, figure 16(a) contains all the data 
that are required for most payload calculations. Appearing as a parameter on figure 

16(a) is the velocity relative to the Earth at the sphere of influence V , which is related 

s 

to the chemical stage burnout velocity Vj^ by equation (7). 

The AV data shown in figures 16(b) and (c) have been included mainly to illustrate 
the sensitivity of AV with respect to specific impulse and initial acceleration. Figure 
16(b) shows that the all-propulsion AV is a very weak function of specific impulse for 
values greater than 3000 seconds. Figure 16(c) indicates that the all-propulsion AV is 
about 30 percent greater than the high acceleration AV. 

Figure 16(c) also indicates that trajectories with two separate power phases were 
found to be optimum over a wide range of a^. These trajectories consisted of an initial 
power phase, followed by a coast phase, followed by another power phase, and terminated 
with another coast phase. As a^ is increased, the initial power phase becomes shorter. 
If a^ is made large enough, the initial power phase vanishes leaving trajectories that 
commence with a coast phase. Further increases in a^ cause the trajectories to ap- 
proach the limitii^ case wherein the power phase is initiated about 75 days (single im- 
pulse transfer time) before the end of the mission. Alternatively, instead of holding 
mission time fixed, mission time could be defined in terms of when a power phase com- 
mences (which would eliminate initial coast phases). 

CONCLUDING REMARKS 

The principal conclusion drawn from this study is that a combination of high-thrust 
and low-thrust stages can be merged into a hybrid system offering distinct advantages 
over either system separately. All- chemical systems are quite unattractive for the 
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close solar probe mission in view of their very small payload ratios. All-electric sys- 
tems are also unattractive from a payload ratio standpoint unless power supplies havii^ 
specific powerplant masses (a) of less than 50 pounds per kilowatt become available. 

But a chemical-electric hybrid system can deliver significant payloads with the relatively 
high values of a characteristic of present and near future power supplies. For in- 
stance, a chemical-electric hybrid system boosted into Earth’s orbit by a Saturn IB and 
employing electron- bombardment ei^ines in conjunction with a 100-pound-per-kilowatt 
solar cell power supply could deliver a 1000- to 2300-pound useful payload to 0. 1 AU - 
depending upon the mission time. 

Consideration of nuclear rocket systems does not alter the essence of the previous 
remarks. In the first place, an operational nuclear rocket is not in the near future. Of 
more importance, however, is the fact that if a nuclear rocket were available, then it 
would be the nuclear- electric hybrid system, which should be compared to the all- 
nuclear and all-electric systems. Such a nuclear- electric hybrid system should offer 
the same advantages over the nuclear and electric systems as the chemical-electric 
hybrid system offers over the chemical and electric systems. Also, because of the 
nuclear rocket’s inherently large size and weight, Saturn V class boosters will be re- 
quired to insert a 0. 1-AU solar probe spacecraft into Earth orbit if reasonable payload 
ratios are to be attained. 

Propulsion system selection is also dependent upon the electric power requirements, 
a factor that tends to favor electric systems. The chemical-electric hybrid propulsion 
system calls for intermediate electric power levels - about 30 kilowatt for a 0. 1-AU 
solar probe boosted into Earth orbit by a Saturn IB. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, August 30, 1965. 
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APPENDIX - SYMBOLS 


^0 

initial acceleration of electric 
2 

stage, ft/sec 


structure mass of electric stage, 
lb 

g 

local acceleration of gravity, 

rt 

^t 

tankage mass of electric stage, lb 


ft/sec^ 

P 

total power delivered to electric 

So 

32. 174 ft/sec^ 


thrustors, kW 

Ic 

specific impulse of chemical 
stage, sec 

^0 

initial circular Earth-orbit 
radius, ft 

le 

specific impulse of electric 

^s 

radius of sphere of influence, ft 


stage, sec 


mission time, days 

K 

hardware fraction of chemical 
stage 


Earth spiral time of electric 
stage, sec 

k 

ratio of tankage to propellant 
mass of electric stage 

AV 

heliocentric characteristic veloc- 
ity, ft/ sec 

Mh 

hardware mass (propellant 
tanks, engines, structure, 
etc. ) of chemical stage, lb 

Vb 

V- 

chemical stage burnout velocity, 
ft/sec 

circular Earth-orbit velocity. 

Mo 

initial mass of vehicle in circular 
Earth orbit, lb 

C 

v_ 

ft/sec 

escape velocity (value of V. 

D 

Mp 

propellant mass of chemical 
stage, lb 

c 

that results in V. = 0), ft/sec 
s 

initial heliocentric velocity. 

nij 

final mass, lb 


ft/sec 

“L 

payload mass, lb 

Vs 

velocity at r^, ft/sec 

"^0 

initial mass of electric stage, lb 

V® 

Earth’s heliocentric velocity. 

mp 

propellant mass of electric 

ft/sec 


stage, lb 

a. 

specific mass of electric power- 

"^PP 

powerplant mass of electric 


plant, Ib/kW 

stage, lb 

y 

correction factor in spiral equa- 

^PS 

propellant mass of Earth- 


tions 

escape spiral, lb 

V 

e 

thrustor efficiency 
terminal heliocentric travel 


angle, deg 


18 


V 


Subscript: 
opt optimum 


V /I g 
c e®o 

cp terminal configuration angle be- 
tween Sun-probe line and 
Sun-Earth line, deg 


Superscript: 

(• ) differentiation with respect to 
time 
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(b) Using electric propulsion. 

Figure 2. - Maximum payload ratio for 0. 1-AU solar probe. Optimization of payload ratio performed with a /77 held 
constant. 



Thrustor efficiency, rj 



Figure 3. - Comparison between electric system and chemical-electric hybrid system for 0, 1-AU solar 
probe. 
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Figure 4. - Assumed variation of thrustor efficiency as a function of specific impulse. 
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Figure 7. - Maximum payload ratio delivered to 0. 1 AU. Optimization includes thrustor efficiency plotted against 
specific impulse function. 
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Figure 8. - Optimum specific impulse for 0. 1-AU solar probe. 
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Figure 9. - Optimum chemical stage burnout velocity for 0. 1-AU solar probe powered by hybrid system. 
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Useful payload delivered to 0. 1 AU, lb 



Figure 11. - System performance in terms of useful payload, electric power, and 
mission time. Specific powerplant mass, 100 pounds per kilowatt; initial mass 
of vehicle in circular Earth orbit, 32 000 pounds. All systems initially boosted 
into 100-nautical-mile Earth orbit by Saturn TB. 
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Useful payload delivered to 0. 1 AU by nuclear system, lb 
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Figure 12. - Effect of perihelion radius for soiar probe. Hybrid system with 
large electron-bombardment thrustors and 100-pound-per-kllowatt power 



Figure 13. - Trajectory diagrams of two hybrid system missions. 
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Figure 15. - Trajectory diagram of absolute and local maximums. Mission time, 260 days; specific impulse of electric stage, 
8987; optimum travel angle; all-propulsion electric system. Thrust pointers are placed at 20-day Intervals. 
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(a) All-propulsion trajectories. Specific impulse of electric 
stage, 10 000 seconds. 
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Specific impulse of electric stage, Ig, sec 
(b) All-propulsion trajectories. 


1.0 1.8 2.6 3.4 4.2 5.0 5.8 6.6x10’^ 

V9o 

(c) Mission time, 182.5 days; specific impulse of electric 
stage, 10 000 seconds; velocity at r^, 0. 

Figure 16. - Heliocentric AV data. 
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^'The aeronautical and space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administratioyi 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof.” 

— National Aeronautics and Space Act of 1958 
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